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We derive an exact Markovian kinetic equation for an oscillator linearly coupled to a heat bath, 
describing quantum Brownian motion. Our work is based on the subdynamics formulation developed 
by Prigogine and collaborators. The space of distribution functions is decomposed into independent 
subspaces that remain invariant under Liouville dynamics. For integrable systems in Poincare's 
sense the invariant subspaces follow the dynamics of uncoupled, renormalized particles. In contrast 
for non-integrable systems, the invariant subspaces follow a dynamics with broken-time symmetry, 
involving generalized functions. This result indicates that irreversibility and stochasticity are exact 
properties of dynamics in generalized function spaces. We comment on the relation between our 
Markovian kinetic equation and the Hu-Paz-Zhang equation. 
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I. INTRODUCTION 



matrix of the oscillator using a path-integral method, 



A well-known model of quantum Brownian motion is 
a harmonic oscillator linearly coupled to a bath of field 
modes. The Hamiltonian is (with h=l) 
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where qi,pi are the positions and momenta of the har- 
monic oscillator and qk,Pk are positions and momenta of 
the field oscillators. Here k are the wave numbers, and A 
is a dimensionless coupling constant.^ 

The Hamiltonian (|l.ip has been considered in numer- 
ous papers (see Dekker's review [l|). Hu, Paz and Zhang 
have obtained an exact equation for the reduced density 
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are dimensionless coordinates, and the time-dependent 
coefficients are defined in Ref. |[2]. 

The terms with r(t) and T{t)h{t) on the right hand 
side of Eq. (|1.2p suggest the existence of damping and 
diffusion processes characteristic of Brownian motion. 
Strictly speaking though, the Hu-Paz-Zhang (HPZ) equa- 
tion (jl.2p is time-reversal invariant and it corresponds to 
a deterministic evolution of wave functions. Indeed, the 
solution of the HPZ equation is equivalent to the solution 
obtained from Schrodinger's equation, i.e.. 



Pr{t) = Trp (e 
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where Trp means trace over the field and 
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is the initial density matrix, diagonalized in a suitable 
basis of wave functions I'i/'a)- The wave functions \ijja) 
form a complete orthonormal basis of the whole system 
of harmonic oscillator and heat bath. Due to the equiva- 
lence (|1.4I) . the HPZ equation describes a time-reversible, 
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deterministic evolution of each wave function. This con- 
trasts with true Brownian motion, described by a Marko- 
vian equation with broken time-symmetry, which corre- 
sponds to a stochastic evohition of each wave function 

The derivation of irreversible Markovian equations 
from dynamics has been a great challenge This is 
related to the apparent incompatibility between the sec- 
ond law of thermodynamics and time-reversible dynam- 
ics. One point of view is that Markovian equations ap- 
pear as an approximation of the dynamical equations. 
This is the so-called Markovian approximation, valid for 
weak coupling between interacting particles, and for time 
scales of the order of the relaxation time to equilibrium 

However, one can take a different point of view, where 
Markovian equations are formally derived from dynam- 
ical equations without any approximation. This is the 
subdynamics formulation developed by Prigogine and col- 
laborators [1, [1, H B S 113, [iH ■ In this approach, essen- 
tial elements are the distinction between integrable and 
non-integrable systems in the sense of Poincare, and the 
use of generalized functions [l2i] . In this paper we derive 
an exact Markovian equation for the quantum Brownian 
oscillator, based on this approach. This equation is valid 
for both weak and strong coupling. As we will show, for 
weak coupling, it agrees with the HPZ equation. 

A few exact results using subdynamics have already 
been obtained [1, [l^. However, to our knowledge, there 
was no derivation of an exact Markovian equation for 
quantum Brownian motion. Previous formulations were 
centered on density operators. Here we focus on the ob- 
servables, that is, on products of creation and annihi- 
lation operators. This allows us to consider arbitrary 
N-particle sectors in a non-perturbative way. 

This paper is organized as follows. In Section |TT] we in- 
troduce our formulation of subdynamics. As in the origi- 
nal formulation, we introduce the projection super oper- 
ators n and 11= 1 — n, which define invariant subspaces 
of the Liouvillc super operator. In Sec. IIIII we define in- 
tegrability and non-intcgrability in Poincare's sense. In 
Sees. lIVlTVl we construct the H projector for the inte- 
grable case and derive recursive relations for this projec- 
tor. Extending these relations we construct H for the the 
non-integrable case in Sec. IVII This leads to our Marko- 
vian equation in Sec. IVIII In Sec. IVIIII we compare this 
equation with the HPZ equation. Concluding remarks 
are presented in Sec. IIXI Additional calculations are 
presented in the Appendices. 



II. SUBDYNAMICS 

In this section we introduce the main ideas of our ap- 
proach. We focus on the quantum Brownian oscillator 
model. 

We consider a one dimensional space. We start with 
the system in a box of size L and impose periodic bound- 



ary conditions. Then in Eq. (jl.l[) we have 

fc = 27r7i/i (2.1) 

with integer n. We are interested in the limit L — > oo, 
where the spectrum of field frequencies tOk becomes con- 
tinuous. We will assume that Ck = C-k and LUk = i^~k- 
This allows us to to restrict A: > keeping only the sym- 
metric part of the operators, i.e., we set qu — q-k- 
We will assume as well that there is no degeneracy in the 
spectrum of ojk for /c > 0. 

It will be convenient to express the Hamiltonian (jl.ip 
in terms of annihilation and creation operators. We ex- 
press the coordinates qi as 
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where a|, are bosonic creation and annihilation oper- 
ators of the particle {i = 1) and field ({i} — {fc}). These 
operators satisfy the usual commutation relations 



for — 1 or k. For the momenta we have 
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Introducing the notation 



the Hamiltonian takes the form [ij, [l^ 
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where -Evac is the vacuum energy. The interaction has 
the following volume dependence. 



(2.7) 



with Vk independent of L. In the limit L oo, the 
sum over discretized field modes turns into an integral 
and the Kronecker delta function turns into a Dirac delta 
function, 



27r 

T ^ 

k 



dk 



L 
2^ 



5k,k' ^5{k-k'). (2., 



Hereafter, whenever we write summations or Kronecker 
deltas, it is understood that we take the limit L oo 
using Eq. (j2.8p . Also from now on, when we take the 
limit L ^ oo we will keep the energy density of the field 
finite. This means that (with (A) = Tr(Ap)) 



(4afe) 



for L — > oo. 



(2.9) 
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This condition is known as the thermodynamic Hmit. 
We consider density operators that have diagonal (S- 
function) singularities in field-mode representation [10, 
[l^ . An example of this class of ensembles is the equilib- 
rium Gibbs distribution. For these density operators, 

E<4«A^'> - (4«fe) - o{Ly (2.10) 

k' 

Diagonal observables are as important as sums of ofF- 
diagonal observables. Due to this property the separation 
of diagonal and off-diagonal observables, which we con- 
sider below, is well defined in the thermodynamic limit. ^ 
The Hamiltonian in Eq. (|2.6p has the form 

H = Ho + XV (2.11) 

where Hq is the unperturbed part describing free mo- 
tion, and V is the interacting part. Corresponding to 
this Hamiltonian we have the Liouville super operator 
(or "Liouvillian" ) 

Lh = [H, ] 

= Lo + XLv (2.12) 

/,From the Liouville equation, 

i^^p{t) = LHp{t), (2.13) 
we obtain the time evolution of averages of observables 

{W{t)) = i:i{Wp{t)) . (2.14) 

We will consider observables depending only on particle 
operators, expandable in monomials, 

M^ = al"a™. (2.15) 
with m, n > integers. Then we have 

{W{t)) = TT(WPp{t)) (2.16) 
where P is a linear projection super operator defined by 

P\a[ 'a,' [[al a^' 
\ fc=o / 

fc 

for TTii, Hi > 0. This projector singles out products of cre- 
ation and annihilation operators with diagonal field oper- 
ators. Every creation operator present in the product 
has to be paired with the annihilation operator ak- 



The projector P commutes with the free Liouvillian 

PLo = LoP. (2.18) 

Introducing the complement projector Q — 1 — P we have 
PQ — QP = 0. Thus under the unperturbed time evo- 
lution (with A = 0), any density operator can be decom- 
posed into two components that evolve independently 

p = Pp + Qp. (2.19) 

Each component remains invariant under the free time- 
evolution. For A = we have 

d 

i-Pp(t) = LoPp{t), 

i^^Qpit) - LoQp{t). (2.20) 

Each component follows its own subdynamics, with 
closed time-evolution. This separation allows us to calcu- 
late {W{t)) knowing only the Pp component of p, without 
the complement Qp component. 

On the other hand, the interacting Liouvillian (with 
A 7^ 0) does not commute with P. We have 

i^Pp^^) = PI^HPit) = PLnPpit) + PLnQpit) . 
at 

(2.21) 

This is no longer a closed equation for the component 
Pp(t). This is the main problem of non-equilibrium sta- 
tistical mechanics. A common approach to deal with 
this problem is to write a hierarchy of equations of the 
BBGKY type Q- Alternatively, one can try to obtain 
closed non-Mar kovian equations (with memory terms), 
such as the Prigogine-Resibois generalized master equa- 
tion [3]. As shown by Hu, Paz and Zhang, for the quan- 
tum Brownian oscillator it is indeed possible to obtain the 
closed non-Markovian equation (|1.2p for the reduced den- 
sity matrix. The non-Markovian character of the equa- 
tion is manifested in the time-dependent coefficients. 

In the subdynamics approach, we introduce a new pro- 
jector n satisfying the following three conditions 

(A) = n, 

(B) niH = LhH 

(C) n = P + AHi + A^Ha + • • • 

where n„ are independent of A. The last condition means 
that limA^o H = P and 11 is analytic at A = 

Using n we can decompose a density operator into two 
components that evolve independently: 

p = Ilp + np (2.22) 



^ As shown in Ref. 4] condition (A) actually follows from condi- 
^ It is possible to avoid summations altogether, and use integrals tions (B) and (C). In this paper we will verify all three conditions 

from the beginning ITJ. The results are the same. are satisfied. 
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where 11 = 1 — 11. Each component obeys a closed equa- 
tion, 



d 

«— 11/9 = Lnlip, 



(2.23) 
(2.24) 



Hereafter we will focus on the 11 component. As we will 
see, this component gives the closed Markovian equation 
describing quantum Brownian motion. The complement 
component 11 gives memory effects associated with dress- 
ing [H 4. 

We focus on the equation 



d 

i—TT[WIlp] = Tr[WLHllp] 



(2.25) 



In the following, we will derive an explicit form of this 
equation. Using the property 



Tt[A ■ Sp] = Tr [{S^A^Yp] 



(2.26) 



where 5' is a super operator and A is an arbitrary oper- 
ator, we have 

i^TT[Wnp] = Tr [{LhH^'W^)'' p] . (2.27) 

To obtain the kinetic equation, we need to calculate the 
quantity LhU'^W'' with W = al"a™. This will be done 
in the following sections. 

We note that (see Eq. ([216]) ) 

Tt[WLhTIp] = Tt[WPLhTIp] . (2.28) 

As shown in Refs. [1, [13]) we have 

PLHlip = enp, (2.29) 

where is a "collision" super operator satisfying the re- 
lation [9, P] = 0. We get an exact, closed Markovian 
equation 



d 

i— PHp = ePIlp 
ot 



(2.30) 



for the component PIlp. This will be verified through 
the direct calculation of Eq. (|2.25p . 

In the construction of 11 we will consider two cases 
discussed next: integrable case and non-integrable case. 
As we will show the Markovian dynamics of Brownian 
motion occurs in the non-integrable case. Our approach 
will be to first obtain 11 for the integrable case, and then 
extend this result to the non-integrable case. 



In a more detailed formulation of subdynamics (see 11]) both El 
and n are further decomposed into a sum of orthogonal projec- 
tors n = n''^' and fi = n'-"^. Each sub-component gives 
a closed Markovian equation. However, sums of these projectors 
can give a non-Markovian equation, as is the case for II. 



III. INTEGRABLE AND NON-INTEGRABLE 
CASES 

In this section we specify what we mean by integrable 
and non-integrable cases. 

For the integrable case, the P and 11 projectors can be 
related by a similitude transformation. 



n = U^^PU 



(3.1) 



where U is a time independent unitary transformation. 
This is the same transformation that puts the Hamil- 
tonian in a diagonal form with no interactions (see Eqs. 
(14. ip , (|4.2p below) . In this way the interacting system can 
be mapped to a non-interacting system through a unitary 
transformation. We call this case "integrable" because 
there exists a one-to-one correspondence between unper- 
turbed and perturbed invariants of motion. Furthermore, 
the perturbed invariants are expandable around A = 0. 
These properties were studied by Poincare in the context 
of celestial mechanics, so when we speak about integra- 
bility, it is in Poincare's sense [2^ . 

In contrast, for the non-integrable case the interactions 
cannot be transformed away through a unitary transfor- 
mation. There is no more a one-to-one correspondence 
between unperturbed and perturbed invariants. The P 
and n projectors are now related by a non-unitary trans- 
formation A, 



n ^ A-^PA. 



(3.2) 



As shown in Refs. [2l|, [24] the transformation A is "star- 
unitary." In this paper we will construct the 11 projec- 
tor directly, without using the A transformation. Let us 
just make a few remarks on this transformation. Rather 
than transforming away the interactions, A takes us from 
the original representation in terms of bare particles to 
a new representation in terms of dressed particles which 
obey stochastic equations breaking time-symmetry. In 
this new representation the effects of noise appear due to 
the non-distributive character of A with respect to mul- 

tiphcation [u, m m, m . 

For the quantum Brownian oscillator we can have both 
integrable and non-integrable cases, depending on the re- 
lation between the frequency of the particle and the fre- 
quencies of the field modes. 

We assume that the field frequencies uJk take the values 



(3.3) 



Here ujq is the lower bound of the spectrum of uJk for 
k = 0. 

The integrable and non-integrable cases correspond, 
respectively, to the following two possibilities p3| : 



(a) uji < ujo 



ib) 



(3.4) 



(3.5) 
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where 



dk 



(3.6) 



The frequency Wc is a threshold frequency for oji, below 
which the oscillator goes from damped oscillation to un- 
damped oscillation. The intermediate case oJo < t^i < t^c, 
gives undamped oscillations as well. In this case the 11 
and P projectors are related through a unitary trans- 
formation, but this transformation is not expandable 
around A = 0. This intermediate case will not be con- 
sidered here. Interesting phenomena associated with this 
case have been considered in Refs. I25l. [2611. 



IV. n IN THE INTEGRABLE CASE 

We consider now the integrable case (a) discussed in 
section [nil In this case the particle cannot resonate with 
the field modes. The Hamitonian (|2.6p can be diago- 
nalized through the unitary super operator U into the 
following form |27j . 



H = ujiA\Ai -I- y^cjfc 



AlAk + E, 



(4.1) 



where u)i is the renormalized frequency of the particle, 
E'vac is the renormalized vacuum energy, and the opera- 
tors A are the renormalized operators of the bare opera- 
tors a. The A operators are given by the unitary trans- 
formation 



A] = U 'a 
A, = 



(4.2) 



We use bars to denote variables in the integrable case. 
The transformed operators satisfy the relations 



LhA\=QiA\, LhAi 
LhAI = ujkAl, LhAu 



-LUlAl, 

-uJkAk ■ 



(4.3) 



As mentioned in the previous section, we call this case in- 
tegrable because this system follows Poincare's criterion 
of integrability. There is a one-to one correspondence be- 
tween the unperturbed invariants of motion aja^ and the 
perturbed invariants The perturbed invariants are 

expandable around A = 0. 

The super operator U^^ may be written in the form 



u au 



(4.4) 



where m is a unitary operator. Thus we have the dis- 
tributive property 



(4.5) 



It follows that the operators A satisfy the commutation 
relations 



[^., A]] = h 



(4.6) 



where i, j = 1, k. 

The transformed operators are given by the linear re- 
lations 



A 



Ai 



-* t 
"ii'^i 



-* t 



-4iai 



k 
k' 



+ }2dlkak, (4.7) 

k 



k 



4fc'«fe'(4.8) 



with the coefficients c and d written in Appendix \X[ 

The renormalized frequency uji is the solution of the 
equation 

G((I'i)-i=0. (4.9) 

satisfying the condition liniA^o = ^^i j where G is the 
Green's function. 



dk 



'iUlUJkX^vl'^ ^ 



(4.10) 



defined here for lo < ljq. Using the commutation relations 
(l46l) we invert (l47l) to obtain 



-i\ = ciiA\ - dl^Ai 



k 



CkiAl-Y^dl^Ak. (4.11) 

k 



We verify now that in the integrable case, the 11 projector 
is obtained through the relation 

IV = W = U-'^PU, (4.12) 

where we use a bar to remind that this corresponds to 
the integrable case. To prove this, we will check that this 
expression satisfies the conditions (A-C) in Sec, [TTl 

Condition (A) is satisfied, since P = itself is a 
projector. Condition (B) means that 



or 



U-^PULh = LhU^^PU 



PULhU-^ = ULhU-^P. 



(4.13) 



(4.14) 



Using the distributive relation (|4.5p together with Eq. 
(|4.ip one can show that ULhU^^ has the same form as 
the unperturbed Liouvillian Lq, which implies Eq. (|4.14[) 
is true 22]. Condition (C) is satisfied, since the super 
operator U reduces to the unit super operator when A 
0, as can be seen in the explicit forms of the coefficients 
in Eq. (|47l) : see Appendix El 

The distributive relation (|4.5p together with Eq. (|4.12p 
also shows that the projector n can be defined through 
the relation 



— ^1 ^1 li"mfc:nfc^fe ^fc 



(4.15) 



since this is equivalent to Eq. (|2.17p . 

Note that 11^ = II. Henceforth we write nt antici- 
pating the extension to the non-integrable case (see Eq. 

(E2Z1)). 
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V. RECURSIVE RELATIONS 

Before going to the non-integrable case, we will derive 
recursive relations for the n projector in the integrable 
case. Subsequently, these will be extended to the non- 
integrable case as a crucial step in the derivation of our 
Markovian equation. The relations are 

nt (arV/) = ntflt. fit (ar"'ar) 



+ (m-l)Xnt(al™ a'/) 



+ nYU''{a\' 



(5.1) 



and 



fit (apa^*) = nt(apari)-ntai 

+ (n-i)xnt(al™4'-2) 



where 



X = 



^Ckidli{Ak,Al} , 



y 

k 



(5.2) 



(5.3) 



(5.4) 



and { } is the anti-commutator bracket. 

In the rest of this section we present the proof of these 
relations. To facilitate our construction, we define two 
new operators, i.e.. 



(see Eq. (liHt ) and 



dI = CkiAl - dl^Ak 



with their Hermitian conjugates Bi and Dk- 
We have 



Thus 



nt(ar<) 
nt 



(5.5) 



(5.6) 



(5.7) 



(5.. 



+ nt 



k k 

BiiBi+Y,Dir-HBi+J2Dk) 

k k 

Y^Di{Bi+Y^Dir-\B,+j2Dk) 



Using the projection property (|4.15p of n, the first term 
of last expression equals 



nt 



fit . fit 



Bm + Y,Diy"'-\B^ + j2Dk) 

k k 

{B\+Y,DlY--\B,+Y,Dkr 



(5.9) 



Now consider the second term in the right hand side of 
()5.8|) . Expanding this in binomial series we have 



nt 
nt 



^5t.(5t+^i5t)™-i(Bi + ^5. 

m-l / \ ' 



1=0 



l'=Q 



(5.10) 



where C[" =_m!/[(m - l)W.]. We have shifted the D 
freely among B since the D operators commute with the 
B operators (see (|4.6p ). 

Due to nt, in order to produce a non-trivial projection 
the term at the beginning of the product must 

pair up with one of the ^j., D^,, with I possible ways 
or one of the J2k" B>k" with V possible ways. Check- 
ing the volume dependence, we can neglect simultaneous 
pairings of d\. with three or more or Z)^, because 
such terms will be of order 0{1/ L) in the thermodynamic 
limit, and therefore can be dropped in comparison with 
other more dominant terms. For example. 



nT 



L k 



Y,\cki\'AfAl^O{l/L), 



(5.11) 



where we have used Icfcil"* ~ |Vfc|^ '-^ 0(1/ L'^) and the 
fact that {A^i^Af.) ^ 0{L'^) in the thermodynamic limit. 
This last relation follows from Eq. (|2.9p . together with 



(5.12) 



(4ife) = (atafc)+0(l/VI). 



which is due to the volume dependence of the interaction 
Vk- With this consideration, after suitable relabeling I = 



7 



fit 



I - 1 and I' = 1' ~1, Eq. (|5?T0)) becomes 

^, fc' / \ fc / 
E^i^'^1 fe+E^ll 

V fc / \ k / 



m-l-2 



(5.13) 



+ nt 



i'=0 



(m - i)xnt 



+ nYU 



where 



and 



X = 11'' 



E^r 



L fc 



(5.14) 



(5.15) 



which give Eqs. (|5.3[) and (|5.4p respectively. 

With dSSl) and ([SlS]) . we have proven jSH]). Relation 
()5.2|) can be proven in a similar way. 



VI. n IN THE NON-INTEGRABLE CASE 

In the non-integrable case (b) discussed in section [1111 
the particle frequency uji can resonate with the frequen- 
cies LOk of the field. There is no self-adjoint perturbed in- 
variant A| corresponding to the unperturbed invariant 
a|ai, which is expandable around A = 0. This could be 
expected, because now the particle is a damped oscillator. 
Damping comes from resonant emission of field modes. 
Due to damping there is no invariant of the form Ai 
[2^ . The Hamiltonian can now be written as [l^, Il4| 



(6.1) 



where Ak , Aj, are renormalized annihilation and creation 
operators of the field and Evac is the renormalized vac- 
uum energy in the non-integrable case. As in scatter- 
ing theory we can choose either "in" or "out" operators 



|la |. Hereafter we will use "out" operators. As we will 
see, from the "out" operators we will obtain damping for 
i > in the Heisenberg picture. The explicit form of the 
operator Ak is 



Ak = Cfcioi + dkial + 



k' 



Ckk'dk' 



'dkk'al 



(6.2) 



with the coefficients given in Appendix |Al These coeffi- 
cients are proportional to Green's function G~^{ujk) where 



G±(z) = 



dk 



AujiUJkX^vl 



(6.3) 



for general complex argument z with Re(z) > luq. The 
+ (— ) sign means the function is analytically continued 
from the upper (lower) sheet of z. The function G^(z) 
has a pole on the "second sheet," obtained by analytic 
continuation form the upper to the lower half plane of z 
across the branch cut on the positive real axis. Denoting 
this pole as 



Zl = UJl — ij 



(6.4) 



(with 7 > 0) we have G+(zi)^^ = 0. This pole reduces 
to UJl when A — > 0. 

By extracting the residue at this pole in Eq. (16. ip we 
obtain the complex spectral representation (see ^15] ) 



H = ZiAlAi + J2 ^kAlAk + -Evac, 



k=0 



where 



Ak = ^fc [1 + 27ri(wfe - zi)Sc{uJk ~ zi)] 



(6.5) 



(6.6) 



Sc is the complex (5-function, and 

ziA\Ai = -Residue(tjfcij'ifc)^,=^, (6.7) 

00 

= - ^ L0kA\Ak2m{ujk - zi)5c{^^k - Zl) 

k=0 

To evaluate the complex-delta function we first go to the 
continuous limit so the summations go to integrals. Then 
we deform the integration path to a small contour sur- 
rounding the pole Zl. 

By separating the residue at the pole zi we obtain 
the particle operators A\, Ai in the complex spectral 
representation of the Hamiltonian. In this way we obtain 
a closer correspondence between the integrable and non- 
integrable cases. Note that the complex-delta function is 
a generalized function. 

In terms of the non-unitary transformation A men- 
tioned in Sec. IIIII we have 



4 = 

A,. 



A- 



- At a 



A" 



A,. 



t 

i ■ 

A^Oi 



(6.8) 
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The explicit forms of the new operators in Eq. (|6.5p 

are 

A\ = ctial +dtiai+^cjfc4+^dtfeafe, (6.9) 

k k 

Ai = cl^ai+d\^a\+^c[kak+^dl^,al, 

k k 

Ak = Ckiai + dkia\ + ^ Ckk'ak' + ^ dkk'a\' ' 

k' k' 

with the coefHcients presented in AppendixlA] The trans- 
formed operators satisfy the relations 



LhA\ 


— 7* it 

— Zi/l^, 




= -ziii , 




LhA\ 




LhAi 


= -zlA^ , 




LhAI 


^^kAl 


LnAk 


= -LOkAk ■ 


(6.10) 



Due to the complex-delta function, the operators in 
Eq. (|6.9p do not preserve the Hilbert space. For example 
one can show that (see Appendix [B|) 



[A,,A\]=0. 



(6.11) 



and similarly — 0. 

Provided the test functions for integration do not con- 
tain singularities at iVk = z\ or Wfe = ^J", the new set of 
operators obey the commutation relations [lB| 



3.12) 



Ak.Al, 



— 5k,k' 



Other commutators are zero. If the test functions con- 
tain singularities, then we need a careful consideration 
[l^ . Two examples are presented in Appendix |B] Here- 
after we assume the density operator p gives no such 
singularities. 

From Eq. (|6.2p we have 



where 



(6.13) 



(6.14) 



Separating the poles at tOk = zi, zl from Eq. (|6.13p we 
get 



where 



Bl^cnAl~dl,Ai, 



Dl = CkiAl - dl.Ak , 



(6.15) 

(6.16) 
(6.17) 



Note that from Eq. (|6.13p we can calculate the exact 
time evolution of a\ as 

g.L^t^t ^ ^ (dkie'^^'Al - dl,e-'^''Ak) . (6.18) 



From Eq. (|6.18p we can calculate the exact time evo- 
lution of any observable associated with the particle, for 
example its energy. Our goal though is to extract the "ki- 
netic" part of the time evolution, which follows a closed, 
exact Markovian dynamics. This is why we introduce the 
projector 11 (or 11^). 

We will calculate the explicit form of the II'I' projector 
acting on products of creation and annihilation opera- 
tors by extending the results of the integrable case. As 
in the integrable case, the projection II'I' (a|'"a") should 
keep terms where the creation operators A^, are paired 
with the destruction operators Ak- At the same time, 11^ 
should leave intact particle operators A\ and Ai. As we 
have seen, the latter are residues at the poles ujk = z*,zi. 

To define II'I' we start by writing (see Eq. (|6.13p ) 



(6.19) 



We decompose this into a sum of all possible pairings. 
For example for m = n = 2 we have 

af a? - E '^lADkA (6.20) 

ki ■ ■ ■^4 

+ E 'DiDlAsK+ E 'DiDiAA + 

ki , ks ,k4 ki ,k2 , ^4 

klM 

where the prime in the summations means that no sum- 
mation variables are equal. Once we have done this sep- 
aration we can extract the poles of the unmatched oper- 
ators using Eq. (|6.15p . For example we have 



E 'DlDlDk.Dk. (6.21) 

kl ,/C2 ,^4 

= Y^{B\+Y.DlMDkAB,+Y.Dk.)- 

k2 kx k4 

To get the 11^ projection, we simply drop the unmatched 
field operators. Thus we have 

E X<^fc2^'^4-E^i^<^'^-^^i-(6-22) 

/Ci,/C2,fc4 k2 

Since the operators with different ki in the left hand side 
of Eq. (|6.2ip commute, the operators with different index 
in the right hand side of Eq. (|6.22p also commute. 
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In general we can write 

nt(at™<) = nt {B\ + ^ DlriB, + ^ i?.)"(6.23) 

A; fe 

where the projection in Eq. (|6.23p is defined as follows 

\ k=Q ) 

= it™^i-n'^™--4""^r: (6.24) 



which corresponds to Eq. (|4.15p in the integrable case. 

In the Heisenberg picture, Eq. ()6.24|) decays for t > 
and mi, n\ > 0. If we had started with the "in" operators 
we would obtain decay for t < 0. For mi = rii = 0, Eq. 
(I6.24p remains invariant. 

^From Eq. (|6.ip we see that H = H, hence 



Tt{Hp) = Tr(HUp) 



(6.25) 



which shows that pcq = Ilpoq for the equilibrium distribu- 
tion. Similar to Eq. (j6.25[) . for any invariant observable 
/ we have 



Tr(/p) = Tr(/np) 



(6.26) 



Therefore the state Hp contains the total energy and 
probability of p. For non-invariant observables, Hp ex- 
tracts the purely exponential terms of the time evolution, 
which are associated with the complex energies zi,z\. 
In this way, Wp gives the Markovian dynamics of ap- 
proach to equilibrium. The complement component Tip 
extracts non-exponential terms that give memory effects. 
This non-Markovian component contains no net energy 
or probability. 

Following the same steps as in the integrable case we 
obtain recursive relations corresponding to Eqs. (|5.ip 
and fO)) , 

+ nrnt(at'"a5'-i), (6.27) 



+ mynt(a{™-ia^) 



where 



X = 



(6.28) 



(6.29) 



(6.30) 



The n super operator we introduced satisfies all the 
conditions of Sec. [Ill The validity of condition (A) is self- 
evident from the identification of 11^ with the projector 
in (jOi)) . 



The second condition (B) can be verified from the re- 
lations 



LhA = zA, li^A^iA, 



3.31) 



with 



k 

z = mizl - nizi + y^(™fc - nk)uJk , 

k 

C = 0,1. (6.32) 

This implies that Tr(y^[L//, n]p) = for all observables 
A and their linear combinations. 

The third condition (C) is verified in Appendix [Cl using 
the recursive relations for 11^. 



VII. EXACT MARKOVIAN KINETIC 
EQUATION 

With all the previous preparations, we are ready now 
to derive the explicit form of the Markovian equation 
()2.27p . As we saw in Section HH in order to obtain this 
equation we need to evaluate Lh^I^'W'' = L//n'''(a|'"a"). 

^■^From the recursive relation (|6.27p we get 

LnUHaY^a^) (7.1) 

+ (to - l)XLHn^{a\"'-\'l) + nYLHn\a\"'~\'l-^) , 

where L^X = = L^Y since X and Y are diagonal 
in the transformed creation-annihilation operators of the 
field. 
Writing 



A\ = ^{c*nBl + dl.Bi), ii = ^(cnSi + dnS^), 



where 



A ^ \cn\^ - \dn\^ = \N\ 



2^1 



(7.2) 



(7.3) 



and using (|6.10p . we find that 



LhBI ^ a*Bl+pBi, 
LhBi = -f3*Bl-aBi, 



(7.4) 



where 



21^1 



«7, 
17 . 



(7.5) 
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/^From Eq. (|7.ip we then infer that 

LHTl\a\^a^l) (7.6) 
= {ma* -na)I{\a\"'a1) 
+ m(in\a\'^-\'l+^) - n/3*nt(a|™+iar') 

- mn [(a* - a)Y + (/3 - P*)X] n\a^;''~\'l~^) 

- m(m- l)[a*X + /3r]nt(al""'a^) 
+ n{n-l)[aX + l3*Y]n\a\"'a'l~^). 

The consistency between Eqs. (I7.ip and (|7.6p can 
be proven by first inserting the recursive relations for 
nt(aj'"a5'), i.e., Eqs. (|07|l and (jOS)) . together with 
Eq. ([7^ into the right hand side of Eq. ([7?T|) . Then we 
verify that the coefficients of II'I' (a|™a") for all values of 
m and n on the right hand side of Eq. (|7.ip are identical 
to the corresponding coefficients in Eq. (|7.6p . 

Now we come to our kinetic equation. We have the 
relations 

nt(a|'"a^X) = X ■IlHa\"'aJl). (7.7) 

because F and X are diagonal in the field operators, and 
is a projection of the diagonal component of the field 
operators. 



Using (|7.7p . we write (|7.6p as 

Lffnt(af"a^) = ntS (7.8) 

where 

^ = (ma — na)(a{ ai) 

+ m/3(aI"-V^+i) - n/3*(4"+Vr') 

- mn [{a* - a)Y + (/3 - /3*)X] (aj^^'a^'-i) 

- m(TO - 1) [a*X + (aj""^a'/) 
+ n{n ~ 1) [aX + /3*r] (al^aT^) • 

Defining 

P = np, (7.9) 

the kinetic equation (|2.27p becomes 

i^Tr[W^p] = Tr[Stp] . (7.10) 

Using the identities for the trace in Appendix [Fl and 
defining ur = Re(a), I3r = Re(/3), we obtain 



i^Tr[W^p] =TrW- laR ([a\a^rp\) 



+ i 



il{Y - X + i) (^[aip, a|] + [ai,pa|] + [a|p, ai] + [a\,pai 

+ /3_R, (^[a|,pa|] - [ai/5, ai]^ +17 (^[a|,pa|] + [aip,ai\j 

+ (aflX + /3_RF) (^[a|/5, a|] + [a|,pal] - [aip, oi] - [ai,pai]^ 

+ i7(F — X) ( [a|p, a\] + [a|, paj] + [aip, ai] + [ai, pa 



(7.11) 



In this form we can already see an interesting property of field components 
this equation: it is closed for the component Pp, as men- 
tioned after Eq. ^T^. To see this, note that Tr[W^p] = 
Tr[W^Fp]. Moreover we have [P,Y] ~ [P,X] - Oil/^fZ). 
Thus in the thermodynamic limit we can move P past X 
and Y . By its definition, we can also move P past ai, a|. 
So in Eq. (|7.11[) we can replace p by Pp. This shows that 
this is a closed equation for Pp. 



The equation takes a simpler form if we assume that 
initially the density matrix is factored into particle and 



Pp{Q) = pM x\{pk{Q) . 

k 

To see this, we write Eq. (|7.11l) as 

i^Tr[Wp] = Tr[W^e'p] 

with the formal solution 

Tr[l^p(t)] = Tr[W^e~^^*p(0)] . 



(7.12) 



(7.13) 



(7.14) 
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The collision super operator has the following operator 
dependence 

9^e{X,Y,a\,ai). (7.15) 

Following an argument similar to the one above Eq. 
(|5.1H) (see also Ref. [ll]), we may neglect correlations 
for Y and X products, i.e., 

= (r)"(x>™ + 0{i/L) (7.16) 

where using Eqs. (|5.12[) and (|7.12p wc have 

{Y)=TrF{Yl[p,m. (7.17) 

k 

and similarly for X. Thus in Eq. I|7.14p we can replace 
X and Y by their initial averages: 

d = 9{(X),{Y),alai). (7.18) 
I 



This means that neglecting 0{l/L) terms, only the par- 
ticle component of the density operator evolves in time, 



pp{t)^pi(t)^l[pkm (7.19) 

k 



and the field density operator drops out after taking the 
trace. Using this result and defining the reduced density 
matrix 



piixi,x[) = {xi\pi\x[) (7.20) 



we write the kinetic equation in the dimensionless coor- 
dinate representation as (see Appendix [G]) . 



. d , 



(7.21) 





f 92 


92 \ 




[dxj 


dx'l) 



+ ijiiY) - (X) + -) 



9 9 



dxi dx'l 



+ i 
+ 



7 
2 

2 
.7 



(xi-a;'i)2 

9 9 



(9^ + 9^) + 



/ I ^■2 I 9 9 



-1- (xi - x'l) 



9 9 



9 9 



+ (aii(X)+/3j^(y)) 



dxi 9.-j + (,9^T + 9^j^"^+"'^^ 
9 9 



2(a;i - x'l) 



+ ^J{{X) (Y)) 



dxi dx'l 
, I I 9 9 



Pi . 



since it is true when averaging over arbitrary polynomial W . Recollecting terms we finally have our exact Markovian 
equation 



9 , 
*9t^^ 

WW 2 ,2,. iOi f 92 92 



(7.22) 



2cji 



-{Xi - X i)p 



2 \dxl dx'l 
9 

9x1 9a;' 1 



Pi 



9 9 

Z7(a;i - x'l) ( — - ) pi - 2i^K{xi - x'lfpi 



uJiJ{xi - a; 1) - — 



9x1 dx'l 



Pi 
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The coefficients are given by 



K ^ {Y)-{X) + -., 
J = Kj'-K, 



J' = {Y) + {X) + -. 



(7.23) 



Using the identities \\A 



k k 

(7.24) 



we get 



K 
J' 
J 



= 4 



Jo 



2, ,2 



(7.25) 



4 / dfcA^z;^|G+(^fe)|X nfe 
Jo 

/"OC 

4/ cifcAV|G+(^fc)P(|zip-a;D 

JO 



where Uk = {a\ak). 

Note that for the integrable case, with no particle-field 
resonance, we have 7 = and the damping and diffusion 
terms in Eq. (|7.22p vanish. Still, the last term remains. 
This corresponds to the excitation of the particle due to 
virtual processes. 

We also note that the kinetic equation (|7.11[) is valid 
for all field distributions, while Eq. (I7.22p is valid for field 
distributions of the form (j7.12p , which are more general 
than Gibbsian distributions. 



VIII. COMPARISON WITH THE 
HU-PAZ-ZHANG KINETIC EQUATION 

In this section we will obtain an expression for the ki- 
netic equation correct up to O(A^), in the weak coupling 
limit A ^ 0. We then compare this to the A^i-limit of 
the HPZ equation and show that they are identical. The 
X^t limit means that we take A ^ 0,t ^ 00 with }?t 
finite. Physically, this means that we consider times of 
the order of the relaxation time t^^x ^ 1/A^ in a weakly 
coupled system. 

We start by assuming the harmonic oscillator is in- 
teracting with a bath of field modes that is in thermal 
equilibrium, and the number density of the field degrees 
of freedom rife satisfies the Planck distribution. We then 
have 

1 11 

'^k + 



^phujk - 1 ^ 

icoth f ^ 
2 V 2 



(8.1) 



where j3 = l/ksT and is the Boltzmann constant. 

In what follows, we are going to use the approximation 
of the function G^{ujk ) shown in Appendix [El First we 
consider the coefficient J, 



f°° 1 
J^-^l dkX'vl\Gt\^iwl-\zx\')in, + -). (8.2) 

^From Eq. ()E9|) we find that in the zeroth order ap- 
proximation 



X\1\G+M\^cjI-cjI)=0 



(8.3) 



The next order can be found using the expression (see 
Eq. (lE5l) ) 



(8.4) 



which gives 



X'vt\G+{uj,r{u:t-u:i) A^ 



Thus we have 
J ~ -2 







dkX^vlV ^ 2 



X'viV — 



coth 



(8.5) 



V 2 



(8.6) 



For the coefficient K we have using (jE9p 

K ~ 2wf / dcjfc 5{ujk - t^OcJfcCoth ( ^^'^ ^ 



icoth f ^ 
2 V 2 



■7) 



Therefore, in the A -approximation, the kinetic equation 
is 



d 



\dxi dx'i^ 
z7C0th(^^^^ {xi-x[)^ 

T^f^^cothf^) 



Now we show that this expression is exactly the same as 
the X^t limit of the HPZ equation. For weak coupling, 
the time dependent coefficients of the HPZ equation ([TT 
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are found to be 

m : 

mhit) c 
mfit) = 



ds r]{s) cos(a;is) , 



1 /•* 

/ rfs ?7(s) sin(a;is) . 

Jo 



ds v{s) cos(a;is) . 



ds i^(s) sin(a;is) 



1 



Jo 



where 



r]{s) 



duj I{uj) sin(a;is) , 



du) 1[LJ) coth(— — ) cos(a;is) , 
2 



(8.9) 
8.10) 

8.11) 

8.12) 

8.13) 

8.14) 
8.15) 



— ~ (^k) 

oo 



k=0 
2-2/ 



(8.16) 



We note that we have added in an extra factor of 1/Mi in 
the expression of compared to due to the different 
definition of the position xi. 

In the X^t hmit we then take t ^ oo. We use the trick 
in [ll[ to evaluate the time integrations. For example we 
have 



c?s sin(wis) sin(ws) = lim£^o+ — / ds 
4 7q 



-1 



1 



i{uji + uj + ie) — i(wi + u — ie) 
1 1 



— UJ + ie) —i{oJi — uj — ie) 
\2i'K5{uJi -\- uj) — 2i7rS{uji — uj)] 

^5{uji - uj) . 



.17) 



The last line is due to the fact that uj can only be positive. 

In this way we can calculate all the coefficients. Using 
the expressions for uji and 7 in Appendix |D] we find 



sn^it) 



dujI(uj)V- 



^ 4 



dkV 



UJlUJkX^vl 



r{t) ~ — / duji{uj)-d{uji~uj) 

= -^I(uJi)=7:X'v\uJi) 
2uj\ 



(8.19) 



V[t)h{t) ~ / dw/(w)coth(/3fiw/2)-^(w- wi) 
Jo 2 

= |/(wi)coth(/3fiwi/2) 

= wi7C0th(/3fiwi/2) , (8.20) 



}? \ dujI{uj)coth{Phuj/2)r- 



10 

~2u)iV 



X2„,2 
dk^^ 



■coi\\{phLJk/2) 



.21) 



Upon comparison of the A^t-approximation of the HPZ 
equation with the A^-approximation of our exact Marko- 
vian kinetic equation (j8.8p . we find that they are identical 



Beyond the weak coupling limit, our equation gives the 
Markovian dynamics of the quantum Brownian oscillator 
valid even for strong coupling, and also valid for any time 
scale. For t ^ 00 the solution of our equation gives the 
equilibrium solution of the complete dynamics. This is 
so because the equilibrium distribution is a function of 
the Hamiltonian, and any function of the Hamiltonian 
belongs to the 11 subspace (see Eq. (|6.25p ). The comple- 
ment component Tip gives all the memory effects, which 
vanish for t 00. 



IX. CONCLUDING REMARKS 

The example presented in this paper shows that ir- 
reversible Markovian dynamics can be regarded as an 
exact dynamics taking place in the subspace of density 
operators Hp, for non-integrable systems in the sense of 
Poincare. The breaking of time-symmetry in the equa- 
tion 

i^np = LH\ip. (9.1) 

at 



~2 1 
UJ-i — UJ^ 



(8.18) 



^ In Ref. [Til the operators P^"^) were used to construct the sub- 
dynamics. As a result, the "anomalous" diffusion term in the 
HPZ equation did not appear in the kinetic equation derived 
in Ref. [Till , because this term belongs to the "non- privileged" 
components. In order to obtain this term one has to include 
"non-privileged" components. In contrast, in the present paper 
we use the operator P = X]^^ p(''i) and the anomalous diffusion 
term is included in the privileged components. A detailed calcu- 
lation shows that both the HPZ equation and our equation are 
consistent with Eq. (136) of Ref. Illl in the one-particle sector. 
The anomalous diffusion term involves higher particle sectors, 
which were not considered in [TTIl. 
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is due to n being non-Hermitian, and appears before 
we take the trace over the field. Thus from our point 
of view of irreversibihty, rather than a consequence of 
coarse graining, is a property of the invariant subspaces 
of the LiouviUian. Time-symmetry breaking appears be- 
cause the construction of 11 involves generalized creation 
and annihilation operators {Ai, Ai). These are eigen- 
operators of the LiouviUian with complex eigenvalues (zi , 
zl) either in the lower or upper half-planes. In this for- 
mulation we add no extra dissipative terms to the Liou- 
viUian. 

The formulation presented in this paper links stochas- 
tic processes and dynamics in a direct way. Once we have 
a Markovian kinetic equation we have a stochastic pro- 
cess described by Langevin-type equations without any 
memory terms. An interesting question is to see what 
is the spectrum of quantum noise associated with such 
Langevin equations (see also [28j). 

In this paper we focused on the lip component of the 
density matrix. In a sense, this component corresponds 
to traditional thermodynamics. From the Markovian ki- 
netic equation we can derive a non-equilibrium entropy 
and the second law of thermodynamics even for strong 
coupling. This could be considered in a subsequent pub- 
lication. 

In contrast to Hp, the complement component Hp gives 
"non-traditional" thermodynamics including memory ef- 
fects. Deviations from thermodynamics in small quan- 
tum systems have been reported in Ref. [2§]. It would 
be interesting to see how the behavior of lip is related to 
these deviations, and what type of non-Markovian equa- 
tion is obtained for this component. 

The model we considered is exactly solvable. For sys- 
tems with non-linear interactions we have to use a per- 
turbative approach. It is our hope that some of the ideas 
presented in this paper will be useful for these systems. 
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APPENDIX A: COEFFICIENTS OF 
TRANSFORMED OPERATORS 

For the integrable case we have 



cn = -TV. 
cifc = N- 



2wi ' 



UJk - UJl 



dn = -N 

dik=-N _ , 

Ckk ^ I , 

Cfel = -AVfcG+(wfc)(wfe +uji) , 
dki = -AVfeG+(wfe)(wfe - LOi) , 



(AI) 
(A2) 
(A3) 

(A4) 

(A5) 
(A6) 
(A7) 



UJk' ~ujk~ie 



{k + k') 



dkk' = -2wiAVfcG+(wfe)- 



XVk' 



UJk' + LOk 

The normalization constant N given by 

AujiUJky^^Vk 



UJl 



1 



dk- 







For the non-integrable case we have [ij 



cii = -N'' 
cik = N 



2uji 
XVk 



du 
dik 

Cfcl 

dki 

Ckk' 



-N 



(uJk ~ zl). 

zl - UJx 



2^1 
\Vk 



UJk + zl 
-\VkG^{ujk){uJk + wi) , 
-AVfeGj(cjfc)(cjfe - cji) , 
\Vk' 

2uji\VkG+{ujk) 



dkk' 



-2uj^\VkG+{ujk) 



UJk' - UJk - le 
XVk' 



Wfe' -I- UJk 



and 



Cfcl — Cfcl, Ckk' — Ckk' 

dki — dki, dkk' — dkk' 



We define Gj(wfc) as 



Gj(wfc) = G^{ujk) - i'^^^^SciuJk 

UJl 



(A8) 
(A9) 



. (AlO) 

(All) 
(A12) 
(A13) 

(A14) 

(A15) 
(A16) 

(fc ^ k') , 

(A17) 

(A18) 



(A19) 
(A20) 



Zl), (A21) 
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where — iV^/(2a;i) is the residue of G~^{u}) at the pole zi 
in the second sheet. The normahzation constant N is 



1 



(A22) 



These coefficients give the "out" eigenoperators of the 
Hamiltonian, A^,, A\ and their Hermitian conjugates. In 
previous pubHcations (e.g. the "in" states were used 
to obtain decaying states for t > in the Schrodinger 
picture. In this paper we consider observables in the 
Heisenberg picture, where the "out" operators are the 
ones that decay for t > 0. 



APPENDIX B: COMMUTATION RELATION 

In this Appendix we prove the commutation relation 

[A,,A\]^o, (Bi) 

Using the explicit forms of Ai , A\ wc obtain 



[A,,Al] = \N\^Y. 



Clfc 



\k=0 



k=0 




(B2) 



We will show that the expression inside parenthesis van- 
ishes, 



We know that zi =0)1 — 17 is the pole of G'+(u;), defined 
in (16.31). Therefore 



~ ^ (ojI~z^) + 

k=Q ^ 1-^ + 



^^^{cok-z^)+ fr'o^k + z^ 

Subtracting the complex conjugate expression and divid- 
ing the result by 2wi{zi ~ z^) we have 



E 



_ zi+ zl _ Cji 

2wi LUl 



LLlk+ Zl\ 



(B5) 



which is equivalent to Eq. (jB3p . This proves the desired 
expression (|B1I) . 

/From Eqs. ^^AJ\ and ([gTS]) we find that [ai,i?j] = 
[^j, Dj:, _bJ]. We then deduce the following other com- 
mutation relations 



Ak.A\ 



2TTiXVk{uJk - zi)5c{uJk - zi) 

G+{LOk) 



-2TriXVk 



(ujk - z*)6c{ujk - z*). 



(B6) 

If the test functions contain singularities at Uk = z\ or 
Wfc = zl, then these commutators are non- vanishing. 



APPENDIX C: PROOF OF ANALYTICITY 

In this section we verify condition (C) on the 11 pro- 
jector for the non-integrable case. This means that 

limnt(al™+ia5') - P(al™+^a5') = at™+^? , (CI) 

A — >-0 

lim nt(al™al'+i) = P(at™a^+i) = al"a^+i , (C2) 

A — >Q 

for all m, n > 0. Furthermore, 11 has to be expandable 
in a power series of A. If 11 satisfies these conditions 
we will say, in short, that it is analytic at A = 0. This 
property is not trivial, because in Eq. (j6.23|) there appear 
non-analytic terms in the products or commutators of 
renormalized operators, as in Eq. (|6.11|) (see also [l^). 
We have, using Eqs. (j6.6p and ()6.15p 



U^a\ = Sj, tfai = Bi. 



(C3) 



Both expressions are analytic at A = 0. Assuming that 
nt(al'"a?), n^al^^-^^) and nt(4"a5'-i) are analytic 
at A = 0, we will show that the recursive expression ()6.27|) 
is also analytic. This will prove Eq. (jCl|) by recursion. 

We start with the first term in the right hand side of 
Eq. (lOTll . 



(C4) 



E^ir 



iB, + J2D,y 



As we show now, this product generates non-analytic 
terms, even if its two factors are analytic. In the prod- 
ucts between outside the brackets and either Bi or B\ 
inside the brackets there appears the term v4|^i, which 
has the following A dependence 

A\Ai — go(A)a|ai + A / dk gi{k, X)a\ak + ■ ■ ■ 



(C5) 



For the perturbation expansion to exist, the functions 
(70(A), gi{k,X),g2{k,X),- • ■ must be analytic at A = 0, 
with (7o(0) = 1. However, 172(^7 A) is not analytic at A = 
0. We have 



A^ / dkg2{k, A)aj.afe 



X^g2{k,X) 
Wl'X'vl 



(C6) 



1 



1 



\N\ 



2JM 

z\ — z 



(Cjfc - Zl)+ (cjfc - z^)_ 

1 



For A ^ the term inside brackets goes to 
1 1 



(a;fc-zi)+ (cjfc-zj')_ 
goes to 

2m6{uk~uJx). (C7) 



Cijfc — [jJi — it cij/c — [jJi + it 
Moreover we have (see Eq. (jD6[) ) 

zx - zl = -2^7 = -2'Ki}?i^(ujx) + 0(}?) (C8) 
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with 



Similarly one can show that 



dk 



This leads to 



dujh 



lim A 52(fc, A) = — ■TT^i'^k - ^i) 
A— >o dk 



(C9) 



(CIO) 



which is non-zero. 

Coming back to Eq. (jC4[) . the term b\ outside brackets 
in the right hand side can pair with either m of the b\ 
or n of the Bi inside brackets. Thus all the non-analytic 
terms involving Jii are included in 



n{B\B,U,,-n\a\^ar^) 

m(Bf)„o„-nt(aI'"-ia5'). 

(Cll) 



In order for Eq. (|6.27p to be analytic, the second and 
third terms in the right hand side of Eq. ()6.27|) should 
cancel the non-analytic terms of Eq. (jClip . Combining 
Eq. ((6?27)l and Eq. ((CTT|) we obtain 



nt(al"+^ar) 



= n(StBi+y)„o„-nt(aI"ari) 

+ m{Bf + X)„on • nt(ai"-V'/) . 

(C12) 



^From Eq. (jClOP we get 



lim(SjBi)„on = - / dujk.5{uJk~ uji)a\ak (C13) 
Jo 

where we used limA„>o |ciip = 1 and limA^o MiiP = 0. 
On the other hand we have 



\cki\'^alak 



dk X^vi\G+{ujk)\'''ii^k +ujifalak 



where we replace ^j^^fc by aj,ak in the thermodynamic 
limit (see Eq. ([5T2)) 1 

In the limit A ^ we obtain (see Eq. (|E9p ) 



lim Yn. 



Thus we get 



1 



dujk —^d{ujk - iOi){u!k + LUifalak 
dwk5(u}k — 0Ji)a\ak ■ (C15) 



\\xa{B\Bi+Y\- 

A^O 



ynon - (C16) 

which shows that the non-analytic terms cancel. 



lim(B|Vx)„o„ = 0. 

A^O 



(C17) 



Since we assume that n'l'(a|™ ^a") and lit (a|'"a" 
are analytic, we conclude that (|6.27p is analytic at A = 0. 
Thus, by recursion nt(a|'"^^a") is analytic for arbitrary 
m, n. 

We can show in the same way that the recursive ex- 
pression (|6.28p is analytic at A = 0, which proves Eq. 
(IC2]). 



APPENDIX D: EVALUATING Coi AND A TO 0{\^) 
^From Eq. (|B4p we have 

(Dl) 



10 i^l - ^?) + 
Approximating ~ o;^ + ie in the denominator we have 

iuJiUJk^'^vl 



Writing 



dk- 



^l-^i-i<^ 



. (D2) 



■'k ^ ^1 



1 c, \ 

= 1^— 2 + d(UJk~UJi), 

we then obtain for the real part and imaginary part of 
(|D2|1 as 

Ql^^l-r r dk ^^^^p^ + 0{\') . (D4) 



(C14) Therefore, 



^k~^i 



Co, ^ uj,-V I dk ^^^^^ . (D5) 
lo ^k - 



and (see Eq. (iU9|l ) 



(D6) 



APPENDIX E: GREEN'S FUNCTION IN THE 
WEAK-COUPLING APPROXIMATION 

In this Appendix we find an approximation of Green's 
function valid for weak coupling. We start by ex- 
panding the inverse function around the pole zi. Since 
this function depends on w^, i.e., 

[G^uDr =^l--l-r dk' t^""f"" , (El) 

Jn ^h' — ^k — 
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we make an expansion in the variable ujf, around z^, 

+ l{u.i~zmGH4rr + ---- (E2) 

We have G^{zf)~^ = 0. The first derivative term is given 
by (with N defined in Eq. ((A22|) ) 



= - 1 



dk 



UJl 1 



(E3) 



For weak coupling we may neglect the second and higher 
derivative terms. This gives 

G^{uJk) = — — ^T^ — 5- + higher derivatives . (E4) 
Lui uji - z{ 



Another useful formula follows from the exact relation 

4^zA2«2(^^.)^^|G+(cc;fc)p - G+K) - G-^) (E7) 

where u^(ti'fc) = v^, dk/du)k. Using (|E5p and (|E6[) . we find 
that 



G+K)-G-(c^fc)~— (SK-u;i). 



Combining this result with Eq. (|E7[) we get 



(E8) 



AX'i\ujk)\G+{Lou)Y 



(E9) 



in the lowest order approximation in A expansion. 
APPENDIX F: TRACE RELATIONS INVOLVING 

al AND ai 



Furthermore for weak coupling we have 1 + 0{X^), 
zi = LOi + 0(A2). Thus we get 



G+K) 



and similarly, 
G-{ujk) 



1 



(Zl - Wfc)(zi +UJk) 

1 

(uJi -Uk -ie){uji +ujk) 
1 



-'P 2^ 2 + iT^Siujk - UJl) (E5) 
uj^ — UJl 2uji 



ujI- ujI+ ie 



Using the following relationships, 



ai,a[ 
'alaT 



a[ai, a[ 
a\ai, a™ 



na 



= —2ma[ 



fn+l m-l 



2nal"~'a5"+i + n(n - l)aT"-"a™ , 



(Fl) 
(F2) 
(F3) 
(F4) 

(F5) 

tn-2 



(F6) 



5{u!k — UJl) . (E6) we can show that 



(n — TO)Tr 
(n + TO)Tr 
nmTr 
Tr 



fn-l m-l ~ 



fn+l m— 1 ~ 



a 



Tr 



Tr 



Jn-l m+l ~ 



I 1 \ t" m— 2 ~ 

raym — lja| p 



Tr l)al""^a5"p 



= Tr 
= Tr 
= Tr 
= -Tr 
= Tr 
= -Tr 

= -Tr 



r t fn mi ~ 

[flj^fli, a{ \p 
r t t" im ~ 



tn m ~ t 



a{ ai alpai 



- 2Tr 
(ra + m + l)Tr 



a{ Qi p 



Tr 



a[ Gi aipal 



alaT]p 



Tr 



al ai (alpal — pa[al) 



r t «1 t™+l ~ 

[a[,aj_]a[ p 



Tr 



a[ a-^ [aipai — aiaip) 



(F7) 
(F8) 
(F9) 
(FIO) 
(Fll) 



2ma[ 



fn+l m-l ~ 



a 



2na| 



fri-l m+l ~ 



Tr 



Tr 



[a\ala\''aT]p 



r t" ml ~ 

[aiOi, a| O]^ Jp 



-Tr 



= -Tr 



a|"a™(2a{pa{ — a{a{p — pa\a\ 



(F12) 



a\"a™{2aipai — aiaip — paiai) 



(F13) 



18 



Furthermore, 



a[aip + pa[ai — aipa[ — a[pai + p = ^ (^^i^iP + palai — 2aipa{ + aiajp + paia\ ~ 2a\pai) . (F14) 



r 



APPENDIX G: COORDINATE 
REPRESENTATION OF AND ai 



Starting from 



(gi - i- 



Pi 



) and Pi = 7^, (Gl) 



for an arbitrary vector \(f)), we find that 
{xi\al 



(a;i|ai|^) 

(0|ai|a;'i) 



7! ('^^ + 9^'^^^''^^ 



J_ / , _d_ 
^/2 r ' ^ ax'i 

V2 V ' dx'i 



(G3) 
x'l) , (G4) 



The ket of the dimcnsionlcss coordinate xi is related to 
Qi by \qi) = (Mio^i)^^'' |a:;i). From the relation Oia| — 
ajai = 1, we also find that 



(a;i|oia||^) 
{xi\a\ai\^) 
(G2) (<^|ai4|x'i) 
{(j)\a\ai\x' i) 



d_ 
dx-j 

_ d 
2 V""'" 9^ 



1 
2 

1 
2 
1 
2 

2 V"^ '" 



2 v"^ ' ^ dx\ 



'\x'i) . (G5) We then deduce that 
I 



X 1 



a; 1 



dx'i 

d 
dx\ 



]x\). 
(G6) 



{xi\[alai,p\\x'i) 





dxj dx'l J 



1 + (Xi - X i) 



p{xi,x[) , 



{xi\{[aip,a\] + [ai,pa\] + [a\p,ai] + [a\,pai])\x'i) = 



p{xi,x'i) , 



(a;i|([aip,al] + [ai,pal] - [alp,ai] - [a[, pai])\x\) 
- (a;i - a; i) - — 



9a; 1 dx'l 



^ + ^ ] (x + X 
dxi dx'l J 



pixi,x[) 



{xi\{[a\p, a\] + [a\,pa\] + [aip,ai] + [ai, pai])\x'i) 



{xi-x\r+ , 



d a 



dxi dx'l 



p{xi,x[) , 



(a;i|([al,pal] + [aip, ai])|a;'i) 



d d \ ( d d \ ( d d 



(G7) 

(G8) 

(G9) 
(GIO) 



(a;i|([al,p4] - [oip, ai])|a;'i) = ^ 



(^?-^1)+(l?3- 



92 92 \ , J d 9 



(Gil) 
(G12) 

p(a;i,a;i) , 

p(xi,a;'i). (G13) 
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